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ABSTRACT 

We present a novel method to investigate cosmic reionization, using joint spectral 
information on high redshift Lyman Alpha Emitters (LAE) and quasars (QSOs). Al- 
though LAEs have been proposed as reionization probes, their use is hampered by 
the fact their Lya line is damped not only by intergalactic H i but also internally by 
dust. Our method allows to overcome such degeneracy. First, we carefully calibrate a 
reionization simulation with QSO absorption line experiments. Then we identify LAEs 
{La ^ 10^^'^ erg and EW> 20 A) in two simulation boxes at z = 5.7 and z — 6.6 and 
we build synthetic images/spectra of a prototypical LAE. The surface brightness maps 
i-^ ■ show the presence of a scattering halo extending up to 150 kpc from the galaxy center. 

For each LAE we then select a small box of (lO/i ^Mpc)'^ around it and derive the 
optical depth r along three viewing axes. At redshift 5.7, we find that the Lya trans- 
missivity T^^^ ~ 0.25, almost independent of the halo mass. This constancy arises 
from the conspiracy of two effects: (i) the intrinsic Lya line width and (ii) the infall 
peculiar velocity. At higher redshift, z = 6.6, where (xhi) = 0.25 the transmissiv- 
ity is instead largely set by the local H i abundance and T^^^ consequently increases 
with halo mass, Mh, from 0.15 to 0.3. Although outflows are present, they are effi- 
ciently pressure-confined by infall in a small region around the LAE; hence they only 
marginally affect transmissivity. Finally, we cast LOS originating from background 
QSOs passing through foreground LAEs at different impact parameters, and compute 
the quasar transmissivity (T^^*-* ). At small impact parameters, d < 1 cMpc, a positive 
correlation between T^^*^ and Mh is found at z = 5.7, which tends to become less pro- 
nounced (i.e. fiatter) at larger distances. Quantitatively, a roughly 10 x increase (from 
5 X 10-3 to 6 X 10-2) of rQS°is observed in the range logM?^ = (10.4 - 11.6). This 
correlation becomes even stronger at z = 6.6. By cross-correlating T^^^ and T^^^ , we 
can obtain a H i density estimate unaffected by dust. At z = 5.7, the cross-correlation 
is relatively weak, whereas at z = 6.6 we find a clear positive correlation. We conclude 
by briefly discussing the perspectives for the application of the method to existing and 
forthcoming data. 

Key words: intergalactic medium - cosmology: theory - diffuse radiation - reioniza- 
tion - numerical simulation 



1 INTRODUCTION 

Lyman Alpha Emitters (LAEs) are galaxies showing promi- 
nent emission in the 2p — >■ Is transition of the hydrogen 
atom resulting in a Lya emission line at Aq — 1215.668 
A. At redshift z ^ 2 the line can be observed in optical/IR 
bands, thus allowing searches of distant galaxies traced by 
such radiation. What is the source of Lya photons? Short- 
lived, massive stars produce large amounts of UV photons 
with energy hi/ > 13.6 eV = 1 Ryd which ionize the sur- 
rounding gas. H-atoms recombine on a short time scale in 
the dense interstellar medium and ~ 2/3 of the recombi- 



nation cascade ends up in Lya photons, powering the ob- 
served luminosity. The rapid production of dust associated 
with star formation process, and most noticeably at hig h 
redshift by supernova explosions (|Todini fc Ferrar"al l200ll '). 
might however lead to important attenuation of the Lya 
emission line. In spite of this drawback, Lya is still consid- 
ered the best available tracer of high-redshift star-forming 
gala xies. This was realized alre ady more than four decades 
ago dPartridge fc Peebledl 19671 ). but a full-scale application 
of the method had to await for technological progresses and 
beca me routinely use d only towards the end of the last cen- 
tury (|Hu et al.lll99i ). 



© 0000 RAS 



2 Baek, Ferrara & Semelin 



Since then, several hundreds of LAEs over red- 
shift z ^ 6 have been dete cted usinp narrowband_Jmag- 
ing a nd spectros copy (e.K. , IHu et al . 1998; Rhoads et^lj 
20031: iMalhotra fc Rhoadd l2004l:~lTaniEUchi et al 2005 



Kashikawa et alj |200b: lOuchi et al.i 
Castellano et al.ll2010l : iPentericci et alJ 



200. 



20111 



l2009l . 



Redshift z = 

6 corresponds the end of the epoch of reionization, the 
last major cosmic phase transition in which the gas turned 
from a neutral to an ionized state. Albeit the general 
features of the reionization process are now understood 
and constrained by available data ([Chou dhurv 



Ferrara 



2006 



Bolton fc Haehneltll2007l : llliev et all 2009: Mit ra et al 



2011bl lal). many physical details as e.g., the initial mass func- 
tion (IMF), escape fraction of Lya, UV, and ionizing contin- 
uum photons, intrinsic source spectral energy distribution, 
and interstellar (ISM) and intergalactic (IGM) medium in- 
homogeneities, remain only vaguely known. There are strong 
hopes that LAEs ca n allow us to clarify these issues as 



sever al authors (e.^.. jMiralda-Escude fc Reeslll998l: ISantoi 
Haiman fc Ceni[2005l:lDiikstra et aLlbOOTI: iDaval et al 



2004 




200S 





20091 ) have pointed out 



However, before we can fully exploit the LAE potential 
in this sense, one has to deal with the fact that reionization 
itself affects the propagation of Lya photons. In fact, in addi- 
tion to the aforementioned presence of dust, interstellar and 
intergalactic H i atoms produce an effective opacity caused 
by their large scattering cross-section to Lya photons. As a 
result, inferring the intrinsic Lya luminosity of a source from 
the observ ed one it often qu ite difficult. Further complica- 
tions arise: [Paval et akl l|201ll ) have shown that a degeneracy 
between the IGM ionization state and the clumping of ISM 
dust is imprinted in the LAE visibility. Stated differently, 
a wide range of IGM H i fractions, xm , can reproduce the 
observed Lya luminosity function with different Lya escape 
fractions, fa, from the galaxy, i.e. increased transmission 
through a more ionized IGM is compensated by higher dust 
absorption inside the galaxy. As a result, disentangling fa 
from the transmissivity through the IGM is very challenging. 

To attack thi s problem we propose here a novel strategy. 
lAdelberger et al.l (2003, 2005) introduced a new method to 
measure IGM metallicity and neutral fraction around a tar- 
get galaxy using background QSOs or galaxies. Such sources 
are separated from the target galaxy by an angular distance 
small enough that their absorption spectra trace the IGM 
around the galaxy. They measured the spatial distribution of 
metals and HI using the absorption spectra of background 
QSOs/galaxies at redshift 2 ^ 2 ^ 3. We exploit this suc- 
cessful method to investigate the transmissivity of LAEs 
through the IGM during reionization. Several ten s of QSOs 
are detected at redshift higher than z ^ 5.7 (e .g Fan et akl 
I2OOII . I2OO3I . l2006al : ISongailalliooi : IWillott et al.ll2007l. I2OI0I) 
and new QSOs are still being discovered ( Mortlock et al.l 
I2OIII ') in the United Kingdom Infrared Telescope (UKIRT) 
Infrared Deep Sky Survey (UKIDSS)Q With hundreds of 
LAEs in the samples a.t z ^ 5.7, it is possible to derive joint 
constrains on the IGM during cosmic reionization. 

The main goal of this paper is to break dust-HI the de- 
generacy affecting the visibility of LAEs by isolating through 
the additional information provided by near-passing QSO 
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lines of sight the IGM contribution to the damping of the 
line. To achieve this goal we will rely on a set of high- 
resolution, radiative transfer cosmological simulations. 

In Section [21 we describe how we fine tune the 
H I density field before running Lya radiative transfer. First 
we obtain the H i density field in a (100 /i~^Mpc)'^ vol- 
ume through UV radiative transfer simulations at redshift 
z = 5.7. We normalize the UV photoionization flux so that 
the synthetic absorption spectra match the observed Lya 
transmissivity deduced from QSO targets. Next, in Section 
[3] we present our LAE model study, a numerical aspect of 
LAE for an individual halo. We present the transmissivities 
for both LAE and QSO and find a relation between them 
at different redshifts (Section [4]). We conclude with a crit- 
ical discussion of the results in Section \S\ In appendix, we 
compare two Lya radiative transfer methods, i.e. a full ra- 
diative transfer and a simplified e^'^" model often used in 
the literature. 



2 NEUTRAL HYDROGEN DENSITY FIELD 

The first step of the method consists in building an accu- 
rate description of the IGM H i density field at the epoch 
at which high-redshift LAE are observed. This is a cru- 
cial step as, due to its resonant radiation transfer of the 
Lya line is very sensitive to the precise value of xm ■ The 
H I distribution is obtained through cosmological simula- 
tions described in the following; however, it is of outmost 
importance that the simulated results are calibrated against 
available absorption line observations before they can be 
used for our purposes. In the remainder of this Section, we 
describe how we take care of this aspect. 

We start from two (modified, for reasons explained be- 
low) snapshots (at z = 5.7, 6.6) of the SI cosmol ogical 
reionization simulation presented in lBaek et al.l (|2010l ) from 
which we take all relevant gas physical properties. Such sim- 
ulation results have been obtained by p ost-processing the 
output of a GADGET-2 (|SpriiLgell [2OO5I ') cosmological hy- 
drodynamical SPH simulation with the radiative transfer 
code LICORICE on an adaptive grid. Such grid has been 
constructed in such a way that each cell contains at most 
TV = 30 SPH particles. This A'^ value gives a minimum cell 
size of 66.88 comoving kpc at z«6, which is smaller than 
the size of a typical, virialized halo at z~6 (~ 200 ckpc). 
The SI run simulates a (100/i~^Mpc)^ volume with 2 x 256^ 
baryon and dark matter particles, which gives a mass res- 
olution of 3.2 X 10^ M© for dark matter and 6.9 x 10* Mq 
for baryons. We assumes WMAP3 cosmologic al parameters: 
Qm = 0.24, ilb = 0.042, h = 0.73, as = 0.76 (|Spergel et all 
l2007h . Reionization starts at z « 14 and ends at z ~ 6 
as shown in Fig. 2 in iBaek et al.l (|2010l '). Dark matter halos 
are identified with a friends-of-friends (FOF) algorithm. For 
each halo, we obtain its dark matter, Mh, and gas, Mg, mass 
and star formation rate ( SFR) M*, with the method de- 
scribed in lBaek et al.|[201ol . All details concerning the source 
modeling (SF R, IMF , SEP , photoionization rate) are de- 
scribed in Baek et al.ll2010l . We keep track of stars formed 
inside each SPH particle according to a Schmidt-law, so to 
at least qualitatively account for photons produced by sub- 
resolution objects. If the star fraction of a particle is >0.1% 
(corresponding to a mass of ~ 10® Mq), we considered it 
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Figure 1. SFR as a function of the halo mass A/^. The red points 
are averaged over Alh, bins span 0.2 dex, represented with the Icr 
error bar. 

as a UV source particle. The number of source particles 
outside resolved halos is about 70% in number, correspond- 
ing to 20% of the total luminosity of ionizing sources. At 
z = 5.7, the simulated mean star formation rate density is 
p* = 8.8 X 10~^ MQyr~^Mpc~^ , in good agreement with cur - 
rent data (|Hopkins fc BeacomI l2006l : iBouwens et al.ll201lD . 
Fig. [T] shows the SFR as a function of the halo mass Mh of 
the galaxies identified as LAEs at redshift z = 5.7. A tight 
correlation between the two quantities is clearly observed, 
with a dependence of the SFR on halo mass approximately 
as oc , with /3 = 1.9 in the range 10.5 < Mh/Mq < 11.5. 
It is assumed that stars form according to a Salpeter Initial 
Mass Function (IMF) in the range 1.6- 12OM0. The galaxy 
Spectral Energy Distribution (SED) has been calculated by 
integr ating the stellar properties give n in lMevnet fc Maeded 
l|2005h and [Hansen fc Kawalei] (|l994 ) over the IMF. We use 
1000 frequency bins between 13.6 eV and 100 eV for stel- 
lar type UV sources, and a photon packet propagates under 
periodic boundary condition. We take into account super- 
nova feedback by simply injecting in the surrounding gas 
particles an amount of 10** erg per unit solar mass of stars 
formed of which 20% (80%) is in thermal (kinetic) form. All 
SPH neighbor particles are kicked with a velocity depending 
on the distance from the center. As we will explain later on, 
matching the observed QSO transmissivity data, requires an 
ionizing photon escape fraction /esc = 0.09. 

2.1 Calibration 

Even a small IGM H i neutral fraction, xhi ~ 10~^ is suffi- 
cient to produce a Lya scattering optical depth 

r(.) = 1.8xlO^(f7M/.^)-/^ (^l±±y\„, (1) 

larger than unity for mean IGM conditions at the redshift 
of main interest in this paper, z = 5.7. Indeed the simulated 
volume-weighted H i neutral fraction is found to be in the 
range —5 < log(a;H7) < —4. Hence, due to the sensitivity of 
T on xm it is necessary to determine the latter quantity with 
a precision < 10"^ The LICORICE code that we use for 
photo-ionizing radiative transfer uses a ray-tracing Monte 



Figure 2. Left panel: Probability distribution function for the 
neutral fraction, xm . Right: Same for the photoionization rate, 
r, computed either directly via full radiative transfer (Fi curve) 
or by assuming photoionization equilibrium for each particle (r2 ) . 
The value of the ionizing photons escape fraction is fesc = 0.09. 

Carlo method on adaptive grid. The snapshot extracted from 
SI shows Monte Carlo noise fluctuations on a;Hi which ex- 
ceed the required precision. To overcome this problem and 
also to normalize the UV flux with QSOs observations by 
varying fesc, we post- process again the SI snapshot. With 
an increased (by a factor 3300) number of photon packets, 
corresponding to 10^" photon packets during a 10 Myr evo- 
lution time; on average ~ lO'' photon packets pass through 
each radiative transfer cell. This procedure enables us to 
achieve a precision on xm and hence to meet the required 
standard as shown in the next paragraph. 

The two panels in Fig. [2] show the probability dis- 
tribution function (PDF) of xm and photoionization rate, 
r, in the high precision run. The xm PDF has a peak at 
Xm ~ 3.2 X 10~^ and closely approximates a log-normal dis- 
tribution. The volume- weighted mean is {a;Hi ) = 2.3 x 10~^. 
At this redshift, the reionization is already completed except 
for very high density clumps (visible as a high xm tail of the 
PDF) where the recombination rate is boosted. From the 
comparison of the F PDFs in the right panel of Fig. [2] we 
can also conclude that the IGM is in photoionization equi- 
librium. The blue curve, Fi, represents the photoionization 
rate distribution for the particles and it has been obtained 
from the full radiative transfer simulation. The red curve 
Fg, is instead the photoionization rate reconstructed from 
the simulated values of xm , gas temperature, T, and den- 
sity, n, assuming ionization-recombination equilibrium, so 
that 

{1 ~ xuifnaBiT) 
l2 = , (2) 

Xm 

where asiT) is the temperature-dependent recombination 
rate. The volume-weighted photoionization rate is (F) = 
1.14 X 10~^^s~^. Both values (xhi) and (F) are consistent 
with the work of iMitra et al.|[2011bl which synthesizes and 
analyzes the observational constraints on reionization. Fig. 
[3] and Fig. [4] show the ihi distributions a,t z = 5.7,6.6, the 
latter obtained by rerunning the snapshot of SI at 2 = 6.6 
with escape fraction f^sc ~ 0.09 for comparison. At z — 5.7, 
the simulation show that reionization is essentially complete 
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and the neutral density field is relatively homogeneous. At 
z = 6.6, instead, large portions (25%) of the simulated vol- 
ume are still essentially neutral and the ionization field is 
patchy with the HII regions confined around sites of intense 
galaxy/star formation (i.e., inside-out topology). 

We caution that the mean xm we obtain at z = 6.6 
might be somewhat high. This is likely caused by the fact 
that in the original SI simulation, appreciable star formation 
appears only around z = 14 (set by the simulation resolu- 
tion). The resulting integrated e.s. optical depth is 0.062, to 
be compared wit h the WMAP7 fid ucial range observed value 
0.088 ± 0.015 (Jarosik et al.llioill '). On the other hand it is 
useful to remember that the Gunn-Peterson effect in high-z 
quasars only gives lower limits for xm ■ The above tension 
does not represent a major issue for our study as in practice 
it only modifies slightly the relation between redshift and 
corresponding mean xm ■ 

The volume-weighted mean {xhi) = 2.3 x 10~^ we 
find is in excellent agreement with the data (e.g. iFan et al] 
l2006bl ) ; however, there is extra information contained in the 
Xm PDF distribution (Fig.[2j left panel) than can be used to 
calibrate the model. To this aim, we use the obse rved Lyg 
absorp tion line spectra of 1 7 z ^ 5.7 QSOs from ISongailal 
112004) an d iFan et al.l (|2006bl ). Let us define the transmission 
at a given redshift, T(2:abs) as the average ratio of observed 
fiux fv to the intrinsic (unabsorbed) one: 

T«^°(zabs) = {.c'vfir'), (3) 

where the average is over a certain wavelength interval along 
the line of sight (LOS). 

In order to compute t'^'^^ from our hydro simulations, 
we interpolate the H i density field on a 512^ fixed grid, us- 
ing the SPH smoothing length of each particle and compute 
the optical depth r. We produce 3000 synthetic spectra by 
piercing through the simulation box with randomly oriented 
LOS, using the same spectral resolution (« 5300) and wave- 
length interval (ISA) to compute Eq. (|3]) as assumed in the 
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Figure 5. PDF of observed/simulated transmission of QSOs at 
z = 5.7 with bin size .02. Filled box plot i s from 17 observed 
QSOs in lSongailall2004l and iFan et"alll2006bl . Dotted box plot is 
from 3000 synthetic spectra of QSOs. 



data (ISongailall2004l ) . The best match to the data is obtained 
by iteration, recursively varying the assumed f^sc value. 

In Fig. O we compare the LOS-averaged simulated 
j'QSO ^j^jj that obtained from 17 observed quasar spec- 
tra. Since the emitted photons are all on the blue side of 
the LyQ resonance, their transmissivity is very sensitive to 
the H I distribution. The form of histogram varies sensi- 
tively with escape fraction; we found that the best fit with 
/esc ~ 0.09 using our source model (initial mass function, 
spectral energy distribution, stellar lifetime etc). We verified 
that for > 100 synthetic spectra t'^'^'-' converges to the fi- 
nal distribution to < 1%. This comparison assures that our 
modeling of the IGM H i density field is accurate enough to 
robustly predict the Lya resonant line transfer from sources 
in the simulation volume. 



© 0000 RAS, MNRAS 000, 000-000 



Joint constraints on cosmic reionization 5 



3 SIMULATING LAES 

Lya photons experience both frequency and spatial diffu- 
sion during the propagation and modeling the Lya spec- 
trum emerging from galaxies requires a high level of ac- 
curacy due to the resonant nature of the Lya line. With 
increasing computational power, several groups have de- 
velop ed Lyg transfer codes (e.glZheng fc Miralda-Escudi 



2002| : lAhn et al.ll2002l:|Diikstra et al.ll2006l: IVerhamme et al 



20061 : iTasitsiomil 12 006': 'Semclin et al.' '2007'; Pierleoni et al 
Laursen et a l. 2009; Fau chcr-Gigucre et al., .201 0). 
Most codes are based on Monte Carlo (MC) methods, ca- 
pable of treating arbitrary geometries and various scatter- 
ing processes. In this Section we describe our LAE model 
with radiative transfer and discuss numerical aspects of the 
LAE modeling. Lya radiative transfer in this work is done 
by LICORICE, whic h has been widely tested and used in 
our previous works llSemelin et all l2007l : iBaek et al] |2009| . 
l2010l ; IVonlanthen et aLlbOllI ). However, differently from the 
ionizing UV radiative transfer, the Lyman line sector of 
LICORICE uses a fixed grid. 



3.1 LAE physical model 

We assume that the main source of Lya emission in LAEs is 
star formation, and compute the intrinsic luminosity of Lya 
from their SFR (see Sec. [2}. About 2/3 of ionizing photons 
are c onverted to Lya photons assuming case-B recombina- 
tion (|Osterbrocklll989l ). so the intrinsic Lya luminosity can 
be expressed as, 



^(1 - fe.3c)QhUa, 



(4) 



where Q is the ionizing photon production rate. The in- 
trinsic luminosity can be attenuated by dust absorption of 
Lya photons in the LAE interstellar medium (ISM). As a 
result, only a fraction /aiL"' can escape into the IGM. Sec- 
ondly, the frequency distribution of Lya photons is broad- 
ened by Doppler effects induced by galaxy rotation. As a 
consequence, the spectrum emerging from the LAE has the 
form 



2 1 

LT[v) = -QhVa{l - fcsc).fa—;=r- 

3 yf-KiXUD 



exp 



(5) 

where /S.vd ~ {vc/c)i'a and Vc is the circular rotation ve- 
locity of the LAE. If star formation is a quiescent process, 
it is Vh ^ Vc ^ 2vh for re alistic halo and disc properties 
||Mo et al.ll998l ; ISantosl2004) . where Vh is the halo rotational 
velocity at radius r2oo within which the collapsed region has 
a mean over density of 200 times the background. We use 
Vc = l-5vh in this work. Finally, due to damping by H i in 
the IGM, only a fraction, T^'*^ of photons escaping out of 
the LAE actually reaches us. Hence, the observed bolometric 
Lya luminosity is. 



(6) 



In the following we will fix fa = 1, i.e. we neglect internal 
dust absorption. This is because our aim is to isolate the 
effects of IGM on the Lya transmissivi ty. A detaile d study 
of the effects of grains can be found in iDaval et al] l|201ll ). 
As already mentioned, they found a dust absorption/IGM 
transmissivity degeneracy that we aim at breaking with this 
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Figure 6. Neutral density field [cm"''] at redshift z = 5.7. The 
comoving pixel size is indicated in each panel. A LAE is at the 
center of the image. 



work through an independent determination of T^^^ . Be- 
fore we proceed, we pause to discuss some interesting side 
issues. 



3.2 Spectral imaging 

In order to produce images and spectra of LAEs, we 
store the data as a thr ee -dime nsional array similarly to 
IZheng fc Miralda-Escudi l|2002l ). Two dimensions corre- 
spond to the sky plane (thus giving imaging information), 
while the third one represents the frequency axis (spectral 
information) like in an IFU data cube. From the beginning 
and at each scattering, we compute the probability, P, for 
the photon to escape along the direction perpendicular to 
the sky plane; the photon is then added to the element of 
the array corresponding to the projected position and fre- 
quency. This procedure is necessary because only a very 
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small fraction of the Monte Carlo photons escape exactly 
in the direction of the image (i.e. observer). The direction 
into which the photon is scattered is deduced from a dipolar 
phase function W{8) oc 1 + cos^O , where 9 is the angle be- 
tween the incident and outgoing directions. The probability 
is added to the array with weight e"^" (1 + )j,^)dft, where dfl 
is the solid angle subtended by the image pixel size, and /i 
is the cosine of the angle between the incident photon and 
the direction perpendicular to the image plane. The optical 
depth is computed from the gas density along the LOS 
to the image, with the frequency v that the photon would 
have if it had been scattered in that direction. The surface 
brightness of each pixel of the constructed image is 



SBr, 



=1 j=i 



k')^le- 



(7) 



7 where L^"' — L'^{v)dv is the bolometric luminosity, 
the luminosity distance to the galaxy, Nph the number of 
photon packets, k the unit vector of the incident photon and 
k' is the one to the observer/image. The sum is performed 
over all photons and all of their scatterings. LICORICE 
uses the comoving frequency and varying expansion factor 
a{t) as a function of the propagation time t of the photon 
iBaek et al.ll2009l ). This method is more accurate than the 
usual approach which uses the same expansion factor in the 
whole cosmological simulation box and avoids first order er- 
rors in {6a /a), which are not negligible when the propagation 
time t reaches several tens of Myr. To determine frequency 
diffusion after each scattering, we take into account pecu- 
liar velocities; the Hubble flow is also included as explained 
above. 

Galaxies with La ^ 10*^'^ergs~^ and observed equiv- 
alent width EW ^ 20 A are identified as LAEs. For each 
identified LAE in the simulation box, we select a smaller 
volume of (10/i~^Mpc)^ centered on it. Such volume size is 
chosen in such a way that all Lya photons on the blue side 
of the line are redshifted and hence not scattered anymore; 
we have checked that on sub- volumes > (10/i"^Mpc)^ T^^^ 
is not seen to vary any longer. We then interpolate all physi- 
cal properties of the gas particles in a grid of 128^ weighting 
on the kernel, w, i.e. a spherically symmetric spline function 
whi ch depends on the smoothing l ength h of SPH simulation 
fsee lMonaghan fc LattanzidI 19851 ). A resolved halo contains 
several tens of SPH particles within 10-100 kpc, so we can 
obtain higher resolution by interpolating on a finer grid as 
shown in Fig. [S] The SPH density and velocity distribution 
can be sensibly improved by interpolating on progressively 
finer grid using SPH smoothing length. We check the to- 
tal transmitted flux varying grid size as shown in Fig. [B] 
Above grid linear size 10/i"^Mpc/128, T^^"^ begins to con- 
verge. The pixel size is 16 kpc (physical) or 2.7" at z = 5.7; 
the probabilities are collected in a spectral array of 240 bins 
spaced by 0.1 A in wavelength, thus yielding a total spectral 
image depth of 24 A (rest frame). 



3.3 A prototypical LAE 

To illustrate the physical features of a prototypical LAEat 
z = 5.7 we select a halo with Mh = 8.6 x 10^" Mq, intrinsic 
luminosity = 7.3 x lO^^ergs"^, and rotation velocity 
Vc = 200 km s~^. These properties are consistent with those 



typically inferred for LAEs. The H i column density inte- 
grated over the line-of-sight through the center of the simu- 
lation box (10/i"^Mpc) is iVm- 6.5 x 10^*^ cm"^. Assuming 
a static gas with homogeneous temperature T — 10* K, this 
Nm value gives a Lya optical depth at line center to ~ 1900. 

Fig. n\ shows an overview of the results. Panel (a) 
is the LyQ surface brightness map obtained by collapsing 
data along the spectral array dimension. In practice, the 
current detection threshold of LAE nar rowband imaging 



is SB « 10 ^*ergs ^cm ^arcsec ^ ( e.g. Shimasaku et al.l 



I2OO6I : lOuchi et alTl2008l : IZheng et af] I2OIGI '). Only the cen- 
tral pixel {SB = 6.6 x 10~^*c.g.s.) exceeds such detection 
limit. This is because we are assuming a point source located 
at the galaxy center and therefore the photons reaching us 
are predominantly those escaping directly with little contri- 
bution from scattered radiation. At larger distances instead, 
all the flux is produced by Lya scattering by H i in the halo; 
its surface brightness dims with radius due to the progres- 
sive decrease of H i density. The halo SB is potentially at 
reach of future instruments and can be used to gain unique 
information about the IGM density structure around the 
emitter. 

The spectra of the source at three different spatial lo- 
cations are also shown in Fig. [T] The black dotted curve 
in each sub-panel shows the intrinsic Lya profile emerging 
from the galaxy before being filtered by the IGM; the solid 
curves refer to the transmitted spectrum. The spectrum of 
the central pixel is very asymmetric as photons in the blue 
side of Lya are completely suppressed as they are redshifted 
through the line core, while photons on the red side are fully 
transmitted beyond A > 8153 A. 

Let us now analyze the spectra in two external regions, 
50 < r < 150 kpc, and r > 150 kpc. The SB observed at 
those locations is due to scattered photons escaping from 
the central pixel. Moving outward from the source, both 
uhi and the magnitude of the peculiar velocity decrease 
with consequent reduction of the frequency diffusion; in ad- 
dition, the damping of the line core becomes less prominent. 
Both spectra show two sharp peaks in analogy with the case 
of an expanding sphere (Somclin ct al. 2007); the left peak 
(blue side of Lya line) is more suppressed than the right 
one because of Hubble expansion, which starts to dominate 
over the infall velocity for distances > 100 kpc. The sepa- 
ration between the two peaks (« 5A) is mainly set by the 
gas velocity field rather than by the H i density since A^'hi is 
relatively small. In fact, with T = 10"* K and tq = 1900, the 
analytic solution of lNeufeldl (|199G| ). holding for a monochro- 
matic source and an extremely optically thick systerr0, gives 
a separation of ~ 1 A at redshift z = 5.7, i.e. smaller than 
what is obtained here. 

Panel (c) of Fig. [7] shows the spherically averaged Lya 
SB profile as a function of galactocentric radius. The SB 
profile steeply {SB oc r"'') drops until r ~ 50 kpc, from 
where it flattens into a plateau extending to r ~ 150 kpc, 
followed by a final descent with slope similar to the inner 
one. The SB dimming is clearly related to the decrease of 
the H I column density with radius. The presence of spheri- 
cally symmetric Lya halos around high redshift sources was 



^ See iDiikstra et al.l l|2006t) for the static, homogeneous sphere 
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Figure 7. (a) hya surface brightness image of a prototypical LAE at z = 5.7. The two circles (of radius r = 50 kpc and r = 150 kpc) 
denote the different regions for which integrated spectra are shown in the central panels. Panels (b)-(d) show the hya spectra, whose 
central frequency is redshifted at Aq = 8150.55 A. The dotted curve in each panel shows the emerging Lya profile, i.e. before any IGM 
damping. Solid curves refer to different regions: (b) central pixel; (c) 50 < r 150 kpc; (d) r > 150kpc. Panel (e) shows the spherically 
averaged Lya surface brightness profile with bin size of 16 kpc. 



pointed out by iLoeb fc Rvbickil (|l999l '). who predicted the 
existence of a roughly uniform surface brightness region of 
size O.lr* = 0.67(ni,/f2A/)pMpc = 100 kpc for the cos- 
mology adopted here. Our simulation, allowing a more pre- 
cise description of the H i density distribution, is in rough 
agreement with the above conclusion if the uniform surface 
brightness is identified with the plateau extension; the SB 
profile is however not constant within the Lya halo as it 
steepens towards the center. 

The above agreement might be somewhat coincidental 
for the following reasons. In our simulation, the radius at 
which the Lya optical depth approaches unity, hence allow- 
ing photons to free-stream away from the source, is r* = 500 
kpc, a fa ctor 2 smaller t h an pr edicted by the analytical 
model of iLoeb fc Rvbickil (|l999l '). which was obtained by 
assuming a completely neutral, average density, homoge- 
neous gas with no peculiar velocity added to the Hubble 
flow. In contrast, our results show that the typical xm in 
the halo is 10~^ — 10~* and the mean gas over density is 
100 — 200. If these values are used, r* becomes much smaller 
than the analytical prediction. However, the further inclu- 
sio n of peculiar velocities, r estores the rough agreement with 
the lLoeb fc Rvbickil (|l999l ) model. 

The previous argument demonstrates that infalling gas 
is the dominant process causing a reduction of transmissiv- 
ity: when infall is artificially turned off, in fact, transmissiv- 
ity ri"*^ for the selected LAE increases from 0.25 to 0.35. 
In principle, in addition to infall, additional peculiar veloc- 
ities might be induced by outfi ows, which can mod i fy the 
visi bility of LAEs. For e xample, iDiikstra fc Wvithd (|2010l ') 
and lDiikstra et al] (|201ll '). by using the spheric ally symmet- 
ric th in shell model introduced by Verham me et all (|2006l . 
I2OO8I ). found that the transmissivity can be ^5 — 10% even 
through a fully neutral IGM. Although outfiows are included 
in our simulations, their effect on visibility appears to be 
rather minor. The only visible effect of the galactic outflow 
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Figure 8. Peculier velocity (relative to the velocity of the central 
pixel) profile along the line passing through the central pixel along 
the line of sight (z-axis). Vertical dotted blue lines represent the 
limit of the hya diffusive area, from which photons enter free- 
streaming region. Velocity profile along z-axis is not symmetric 
because of other halos reside near the l.o.s.. Positions of two halos 
(front/back with respect to the central halo) are indicated with 
blue circles, whose radius show relativistic mass size. 



emerging from our prototypical LAE is to decelerate the 
infalling gas, as it can be seen from the change in the ve- 
locity profile slope within the central 50 kpc (see Fig. |8]). 
Stated differently, it appears that the outflow expansion is 
severely quenched by the ram pressure of the infalling gas, 
thus affecting only mildly the flnal transmissivity. A simple 
estimate reinforces this conclusion. The radius Re at which 
the infalling gas kinetic input rate becomes equal to the 
mechanical luminosity, , of the outflow can be computed 
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from the following equality 

^jL/ = 'inR^.Api„f{z)vfnf = MtvEoe^ = Lm, (8) 



7?e = 70 A" 



-1/2 



Meyr- 



1/2 



"in/ 



lOOkms' 



-3/2 



kpc. (9) 



In the previous equation A is the infalling gas over density 
with respect to the background density Pi„/(z) = A{p(z)); 
we have assumed 1 SN, of total energy _Eo = 10^^ erg, ev- 
ery 100 solar masses of stars formed and a kinetic energy 
conversion efficiency =0.1. It is clear that the infall can 
quench the outflow on very small scales, and therefore its 
influence on the escape of Lya photons is marginal. Very 
likely though the actual situation could be far more complex 
than it is possible to mode here: hydrodynamical instabili- 
ties might perturb the wind/infall interface; the wind might 
be asymmetric and or bipolar. These complications require 
a dedicated investigation that is beyond the scope of this 
study and we leave them for future work. 

The results shown in this Section are based on full radia- 
tive transfer of Lya photons. However, this procedure results 
in relatively large shot noise due to Monte Carlo sampling, 
making the comparison with the QSO transmissivity (see 
below) challenging. For this reason in the following we will 
use the e~'^" model, most often used in literature, to com- 
pute T^^^ rather than the full radiative transfer. We have 
performed a careful comparison (discussed in the Appendix) 
between the two methods, from which we conclude that the 
differences are not significant. Therefore we use the expo- 
nential model as the fiducial case in the rest of the paper. 



4 LYMAN ALPHA TRANSMISSIVITY 

The main aim of this paper is to put joint constraint on 
the LAE Lya transmissivity using information coming both 
from the LAE itself (i.e. in emission) and from the absorp- 
tion spectra of background quasars. To perform this task, 
we cast thousands of LOS around each LAE halo to produce 
synthetic absorption QSO spectra with varying impact pa- 
rameter. By analyzing these two sets of data in combination, 
we show in the following that it is possible to gain key in- 
formation on the environment in which LAEs are embedded 
and to break the dust absorption/HI scattering degeneracy 
that plagues the use of LAEs as reionization probes. We dis- 
cuss the properties of the Lya transmissivity in LAEs and 
QSO separately first; next we combine the two result sets. 



4.1 Lyman Alpha Emitters 

Using the friends-of-friends (FOF) algorithm, we identified 
3070 halos at redshift z = 5.7, corresponding to a halo num- 
ber density n = 1.16 x 10~^Mpc~^. For each halo, as ex- 
plained above, we select a sub-volume of (10/i~^Mpc)"^ cen- 
tered on it and interpolate all physical properties in a grid of 
128^ cells. We use the e model (see Appendix A) and we 
average over 3 LOS along the x, y, z axis from the halo center 
to the edge of simulation box. We identified ~ 2900 objects 
meeting the LAE selection criterion (L^"' ^ lO^^'^erg s"^ 
and EW^ 20A). Fig. [9] shows T^^^ as a function of the 
identified LAE host halo mass Alh- The various radiative 
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Table 1. Radiative transfer parameters. Model are represented 
by curves of different colors in Fig. |9] PV is the peculiar velocity 
of the gas. Vc is the velocity of galactic rotation which determines 
hya line broadening, (see EqO 
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Figure 9. T^^^ as a function of halo mass Mh for identified 
LAEs. The halo mass bins span 0.2 dex and shaded area represent 
the la error bars in each mass bin. Radiative transfer parameters 
for each model is described in Tab[T] 



transfer parameters for each curve/model are described in 
Tab. [U T^*^ depends on several physical conditions, such 
as the LAE luminosity, clustering, Lya profile width, gas 
infall, surrounding whi. The first three properties are posi- 
tively correlated with T^^^ while the remaining ones tend to 
suppress T^*^ ; the net effect therefore involves a complex 
interplay among them. We run several simulations to under- 
stand how these elements affect T^*^ . L3, implementing the 
e"^" model, is our fiducial run; L2 uses instead full radiative 
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transfer. As discussed in Appendix A, transmissivities of L2 
are slightly larger than L3 but show similar pattern. 

Except for the first three low-mass bins, T^^^ — 0.2 — 
0.3, almost independent on halo mass. This is an indication 
of a balance between transmissivity-enhancing and suppress- 
ing factors. For example, massive galaxies have larger SFR 
(see Fig. [T| and consequently higher photoionization rates. 
As a result, their H ii regions are larger and more trans- 
parent to Lya photons (larger T^^^). Ifowever, massive 
galaxies have also more prominent infall rates which tend 
to suppress transmissivity. This physical interpretation is 
supported by the results from run LI, in which we turned 
off peculiar velocity effects. Consistently with the above, 
Tq increases with halo mass. However, we find that the 
positive correlation between T^*^ and Mh is not as much 
driven by the increased photoionization rate, but rather by 
the broader intrinsic line widths. In fact, at z = 5.7 where 
the neutral HI density is low (xhi) ~ f 0~^ and very homoge- 
neous, the ionizing flux from the LAE itself does not alter the 
ionization state of the IGM, and therefore the dependence 
on the halo mass is cancelled. Instead, if we run the L4 model 
without peculiar velocities and constant Lya line width, we 
recover again a flat slope. This is because broader Lya line 
widths allow the escape of more red-wing photons. We con- 
clude that tI^^^ remains approximately constant with halo 
mass because peculiar velocities and Lya line width effects 
roughly balance each other. 

For comparison, we run L5 at z = 6.6 with same addi- 
tional parameters as L3. At this redshift, reionization is still 
relatively patchy and the effects of the HII regions carved 
by the LAEs themselves becomes dominant. Hence, more 
massive galaxies are embedded in larger HII bubbles al- 
lowing a higher fraction of Lya photons to reach us, i.e. 
T}^^^ increases with halo mass. 



4.2 Quasars 

The next piece of information that we need to compute is the 
transmissivity, T^s° , of QSOs whose LOS passes through 
the environment of a foreground LAE at a given impact 
parameter, d (see Fig. llOp . T^^'^is computed following eq. 
[31 but we raised the spectral resolution to R = 40 000 to 
capture properties of the IGM around a LAE as much as 
possible. At J? = 40 000 one spectral bin size corresponds 
to 0.0304 A in the rest frame; for comparison, at z — 5.7, I 
comoving Mpc corresponds to AA = 0.544 A (rest frame) i.e. 
17 spectral bins. We select the transmitted flux of 17 spectral 
bins in the range {Xa ± AA/2)(1 -I- zlae) where zlae is the 
redshift of the LAE. This range corresponds to positions 
in ±0.5 cMpc around the LAE center along the LOS. For 
each LAE, we cast 12,000 LOS along (a;, y, z) varying d, 
and compute averaged values of r^^° in the following four 
impact parameter ranges: d/cMpc G (0 — 0.5,0.5 — 1, 1 — 
2,2-3). 

Differently from T^^^ , at z = 5.7 r^^° shows a strong 
dependence on Mh fFig llip . This effect is caused by the gas 
infall around these halos, whose velocity increases rapidly 
with halo mass, and particularly beyond \og{Mh/MQ) « 
10.8. In addition, r^^° for massive halos {Mh ^ IO^Mq) is 
less dependent on impact parameter than for smaller halos. 

At 2: = 6.6, however, transmissivity drops sharply to 
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Figure 10. Schematic view of a QSO absorption spectrum whose 
LOS passes nearby an LAE. The red (black) spectrum is with 
(without) peculiar velocities. Due to infall, the core of the Voigt 
profile is shifted both to the blue and red side of \a{l + zi^ae) 
causmg icT to mcrease. 
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Figure 11. T^^'^ as a function of halo mass Mfi for identified 
LAEs. The halo bins span 0.2 dex and curves represent median 
of each bin. 
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yet completed ({a;Hi) ~ 0.25) and many smaller halos are 
still embedded in the dense neutral IGM. 

As we have emphasized several times, gas peculiar ve- 
locities are an important physical factor in the determina- 
tion of T^^^ . We define Vinf , the modulus of the infall ve- 



locity, as 
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(10) 
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11.2, since reionization is not 



where Nlos is the number of LOS, Npix is the number 
of pixels around a LAE, Vj is the peculiar velocity of the 
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Figure 12. Gas infall velocity Vi^f in halos of mass M^, repre- 
sented with the Icr (68%) distribution at each bin. 



j pixel and Vi^c is velocity of the central pixel which con- 
tains the halo to the direction of the i — th LOS. We take 
Nlos ~ 3 along the three axis (a;, y, z) and Npix ~ 40 
which corresponds to 3/i~^Mpc (comoving) around the halo, 
i.e. a few times the virial radius of typical LAEs, where gas 
infall should become evident. Fig. 1121 shows infall velocities 
of LAEs as a function of halo mass Mh at redshift 5.7 and 
6.6. Vinf increases with halo mass at both redshifts. Indeed, 
as we will see shortly, «i„/ is a key parameter to understand 
T^so as well as T^^^ . 

Photons from the QSO are in the blue side of Lyo? line 
and they are completely scattered out of the LOS by even 
tiny H i amounts. Usually, the A''hi around a LAE is large 
enough to strongly suppress the quasar flux; however, pe- 
culiar velocities (infall and outflows) can shift the core of 
the Voigt profile and prevent the suppression. Different from 
T^^^ , a large infall velocity increases T^^'^ since the core of 
the Voigt profile is shifted from the central Lya frequency. 
As we show in Fig. IIOI photons in the red (blue) side of 
Aa(l + zlae) interact with the gas having a velocity com- 
ponent parallel (antiparallel) to their propagation direction. 
In both cases, the shift of the Voigt profile core produces a 
T^^° increase. 

Panel (a) of Fig. [T3 shows all data points for r^^° at 
different impact parameters d, while panel (b) shows the 
median of the distribution, both as a function of Vi^f . Val- 
ues are not presented since they overshoot 
the observational precision (10-*). In this case a clear corre- 
lation between r^^° and Wi„/ exists. The slope is sharpest 
at the smallest impact parameter and at large infall veloc- 
ities. This is expected from the simple argument that the 
environmental properties of the gas deeper into the LAE 
gravitational potential, where the component of the pecu- 
liar velocity along the LOS and A^hi are higher (see Figl6]), 
deviate more from the mean IGM ones. At larger impact pa- 
rameters the slope of the median becomes flatter and closer 
to the mean IGM value ~ 0.04 that we obtained in the Sec. 
[2l The medians become independent of d beyond Vi„f ^200 
kms"^, exceeding the mean IGM transmissivity by a factor 
of about 3. 

For QSOs, the peculiar velocity is the most dominant 
factor to determine their transmissivity. r^^° is computed 



in the wavelength range AA — 0.544 A (rest frame, cov- 
ers the 1 cMpc bracketing the center) corresponding to a 
Doppler shift velocity of 135 kms"^. The magnitude of Vinf 
for massive halo is about 150-200 kms"^, well above the ve- 
locity required to shift the core of the Voigt profile out of the 
above wavelength interval. Thus even if the A^hi is highest 
within the smaller impact parameter (0 < d < 0.5 cMpc), 
the infall velocity is large enough to produce an increase of 
r^^° . Since the infall velocity is correlated with halo mass, 
we observe a similar tendency - the independence of T^^'~' on 
impact parameter - in panels (b) of Fig. [TT] and Fig. 1131 

Panel (c) of Fig. 1131 shows similar data for z = 6.6 for 
the closest impact parameter, which shows a steeper positive 
correlation between r^^° and «i„/. At this redshift, 25% of 
the volume is still neutral and some of the halos are embed- 
ded into completely neutral patches. As mentioned, at higher 
redshift the global ionization state of the gas dominates the 
transmissivity over peculiar velocity effects and the size of 
H II bubble plays a more important role on T^^^ than Vi„f . 
The data of r<?^° at z = 6.6 is more scattered than the 
data at z — 5.7 due to the inhomogeneous A''hi distribution 
and the smaller number of resolved halos. The median of 
r^^° below Vinf = 120km s^^ drops very sharply and even 
the median is under 10~* for Vi„f 5^ 75km s"^. We inter- 
pret this trend as due to the fact that halos with small Vinf 
host less massive/luminous galaxies, possibly embedded in 
smaller HII regions leading to a complete r^^° suppression. 



4.3 Lyman Alpha Emitters/Quasars correlation 

In the previous Sections, we have separately analyzed the 
properties of T^^^ and T^^*"* depending on the relevant 
physical parameters. We now ask the question: is there a 
relationship between the two transmissivities ? To clarify 
this point, we plot the r^^° - T^^^ correlation in Fig. [H A 
general property of the two transmissivities is that T^^'"' is 
sensitive to various physical parameters (e.g.infall velocity, 
clustering, neutral fraction) while T^^^is not. This is be- 
cause the transmitted photons from QSOs are all in the blue 
side of the Lya line, hence they experience strong scatter- 
ing in the Lya profile core as they are redshifted by cosmic 
expansion. On the other hand, almost all transmitted pho- 
tons from LAE are already in the red side of the line having 
small cross section to H i atoms when they escape from the 
galaxy. Only photons near to the core of the Lya line are 
affected by infall or H i fraction. 

At 2; = 5.7, we do not find a strong correlation between 
the two transmissivities. Most of the T^^^ data points loiter 
in 0.2-0.3, whereas r^^° shows a large scatter in 10"^ — 1. 
r^^^ is weakly anti-correlated with Vi^f and this erases the 
possible correlation with r^^° , which instead is strongly 
correlated with the infall velocity. However, a,t z — 6.6 there 
is a clear positive correlation between the two transmissiv- 
ities, although the dispersion around it is relatively large, 
as shown in panel (b) of Fig. 1141 Almost rP^° data below 
r^^^ < 0.09 drop under 10"" which is the detection crite- 
ria, but they increase very sharply from 10"* to 0.1 between 
r,^'^'^ ^ 0.2 and T^^^ ^ 0.4. LAEs with low T^^^ are em- 
bedded in the almost neutral IGM patches, and the photons 
from background QSO passing in their vicinity have a very 
high probability to be scattered out of the LOS. 
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Figure 13. T^^*^ as a function of Vi^j . (a) all data points, (b) 
the median of the distribution, (c) all data points and the median 
at redshift z = 6.6. 



5 SUMMARY AND DISCUSSION 

The evolution of the luminosity function (LF) of high-z 
galaxies and, in particular, of LAEs is one of the most 
promising tools to study the interplay between early galaxy 
formation and cosmic reionization. However, before solid 
conclusions can be drawn from such a method, a number 
of delicate aspects must be worked out quantitatively. One 
of the most challenging ones is the presence of a degeneracy 
between Lya line da mping by du s t in t he source and/or by 
intergalactic H i . As^ Daval et al.l (1201 ih have pointed out, a 
wide range of H i fractions is compatible with the observed 
Lya luminosity function since the effects of a largely neu- 
tral IGM are compensated by a higher escape fraction of 
Lya photons from galaxies containing only small amounts 
of dust. With the aim of breaking this dust/H i degeneracy 
we have proposed a possible solution which relies on an in- 
dependent determination of xm exploiting the presence of a 
background luminous source (typically a QSO but the use 
of Gamma Ray Bursts can be also conceived) whose LOS 
pierces through and samples the target LAE surrounding 
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Figure 14. Ta^'~' (all with impact parameter < d < 0.5 cMpc) 
as a function of T^*^ . T^^^ bins span 0.05 and shaded area 
represent la error bars in each T^^^ bin. All T^^*^ data below 
10~* is replace by lO"** which is the actual detectable criteria 



matter. The theoretical basis of the method, which we have 
applied at two different redshifts z = 5.7 and 6.6 where large 
LAE and QSO samples are available, can be summarized as 
follows. 

First we normalize the UV ionizing flux and the H i field 
by using results from QSO absorption line experiments. We 
produce 3000 synthetic QSO spectra at redshift z = 5.7 with 
randomly selected position and direction in a (100/i~'^Mpc)'^ 
simulation volume. Next we obtain the transmissivity for 
each LOS, using the same spectral resolution a nd wave- 
length range as for a s ample of 17 observed QSOs (|Songailal 
12004 iFan et al.ll2006al '). We find that the transmissivity dis- 
tribution changes very sensitively with the mean UV pho- 
toionization rate; the best fit to the data is obtained for a 
volume-averaged (xhi) = 2.3 x 10~^. 

Before running Lya radiative transfer for all identified 
LAEs in the simulation, we have presented the individual 
case of a prototypical LAE to illustrate some relevant physi- 
cal features behind the Lya line transfer. The observed sur- 
face brightness {SB oc r~^) is the sum of the central point 
source and a scattering halo extending up to 150 kpc from 
the galaxy center. Such halo is just below current detec- 
tion threshold and at reach of future experiments. Our re- 
sults also indicate that Lya photons in the red wing of the 
line are largely transmitted at this redshift. We also find 
that infalling gas is the dominant process causing a reduc- 
tion of transmissivity, whereas outflows seem to be efficiently 
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quenched by infall ram pressure, thus confining their efi'ects 
in a relatively small region around the galaxy (< 50 kpc). 

For each resolved LAE in the simulation, we select a 
small box of (10/i~^Mpc)^ around it and derive the optical 
depth r along the three LOS parallel to the three viewing 
axes. At redshift 5.7, we find that Tj*'^« 0.25, almost in- 
dependent of the halo mass. This constancy arises from the 
conspiracy of two different physical effects: (i) the intrin- 
sic Lyof line width and (ii) the infall peculiar velocity. The 
Lya width depends on the galaxy rotation velocity, which 
increases with halo mass; the broader is the line, the larger 
is the fraction of red wing photons that can be transmitted, 
thus increasing T^^^ . The infall velocity also increases with 
halo mass, but it acts to suppress T^^^ , as it blue-shifts the 
line in the rest frame of the infalling H i atoms. Evidently, 
the two effects balance each other almost perfectly, yielding 
a roughly constant T^^^ . At higher redshift, z — 6.6, where 
{xhi) ~ 0.25 the transmissivity is instead largely set by 
the local H i abundance and by the ability of the galaxy to 
carve a sufficiently large H ii region around itself. As the 
SFR and the photoionization rate increases toward large 
masses, consequently increases with halo mass from 

0.15 to 0.3. 

As a next step, we cast thousands of LOS originating for 
background QSOs passing through foreground LAEs at dif- 
ferent impact parameters. Differently from the case of LAE, 
photons from the quasars are in the blue side of the Lya 
line and are scattered out of the LOS by H i in the vicin- 
ity of the LAE. Again, we have emphasized the importance 
of gas infall motions which increase the quasar transmissiv- 
ity (T^^°). At smaller impact parameters, d < 1 cMpc, 
a positive correlation between T^^*"* and Mh is found at 
z = 5.7, which tends to become less pronounced (i.e. fiatter) 
at larger distances. Quantitatively, a roughly 10 x increase 
(from 5 X 10"^ to 6 x 10"^) of T^^° is observed in the range 
logM/i = (10.4 — 10.6). The correlation is even stronger at 
z = 6.6. 

At z = 5.7, we do find a relatively weak correla- 
tion betweenT^ and r^^° . This is because most of the 
T^^'^data points loiter in 0.2-0.3, whereas T^^° shows a 
large scatter in 10~* — 1. However, at z = 6.6 there is a 
clear positive correlation between the two transmissivities, 
although the dispersion around it is relatively large. The 
median is below 10"'' for T^^"^ ^ 0.09, and it mcreases 
very sharply from 10"* to 0.1 between T^^^ ^ 0.2 and 
rpLAE^ 0.4, showing a very strong sensitivity of T^^^to 
T^*^ . The correlation signal is diluted by the fact that the 
properties of the environment of many LAEs (i.e the most 
clustered ones) is not completely determined by the radia- 
tion and peculiar motions caused by the central object alone. 
Thus a QSO line of sight through such a clustered environ- 
ment results in contaminating high T^^° data. Properly se- 
lecting relatively isolated LAEs would boost the correlation 
signal. 

The proposed method therefore appears to be promis- 
ing to determine the evolution of the physical state of the 
IGM with redshift and to allow a detailed study of reion- 
ization. Obviously, its feasibility relies on the availability 
of a sufficient number of high redshift LAEs and QSOs. 
In this sense, we live in an age full of promises. Cur- 
rently, about 60 QSOs a t 5.7 < z < 6.5 have been de- 
tected in total. Of these, iFan et al.l (|2006al ) discovered 19 



QSOs with z ^ 6 and ISongailal l|2004h discovered 5 QSOs 
at 2 > 5.7 in the Sloan Digital Sky Survej|f| (SDSS) us- 
ing moderate resolution spe ctra, R ~ 5000. Additional 19 
QSOs are discovered by Willott et al.1 (|2007| . |2009| . |2010| ) in 
the Canada-France High-z Quasar Survey (CFHQS) with 
moderat e resolution R 1000 - 5 000. In the VLA FIRST 
survejQ, iBecker et all (|2006l . bOld ) have obtained 17 QSOs 
with high {R ^ 40 00) and modera t e (R KibOOO) spectral 
resolution. Recently iMortlock et al.l (|201ll ) discovered the 
highest QSO at z = 7.085. New QSOs are still being dis- 
covered and identifying QSOs at even higher redshifts is 
the next challenge. Ongo i ng ne ar-IR sky surveys like the 
UKIDSS (|Lawrence et al.l [ioOTh will be completed in few 
years. The Large Area Survey (4000 deg^) aims at finding 
z = 7 quasars. Fi nally, the ESO VISTA Telescope surveys 
l|Sutherlandll2009l ) will be even more efficient building on a 
larger camera probably capable of finding z ^ 8 sources. 

On the LAE side, a large sample already exists at 
z 5.7. Up to now 54 (45) LAEs are confirmed by spec- 
troscopic obser vation at redshift z = 5.7 (6.5) in the Sub- 
aru Deep Field (iKashik awa et al.ll201ll ). The number of con- 
firmed LAEs are 70% (81%) of their photometric candidates 
at 2 = 5.7 (6.5). The Subaru/XMM- Newton Deep Survey 
(SXDS) has also found a larg e num ber of LAEs at z = 5.7. 
In a 1 deg^ sky lOuchi et all (120081 ) found 401 photomet- 
ric candidates, 17 of which spectroscopically confirmed. The 
survey locations of SDF (a = 13h24, 5 = 27°) and SXDS 
(a = 2hl8, S — —5°) are within the survey area of VLA 
FIRST and SDSS which provides a first direct opportunity 
to test our method. 

The survey field for LAE is about 1 deg^ in which hun- 
dreds of candidates and tens of confirmed LAE have been 
found. On the other hand QSOs are found over wide area 
over thousands of deg^ and the probability that we can find 
them next to the target LAE is low. A much more efficient 
strategy to find QSO-LAE pairs would be to perform a deep 
spectroscopic sky survey around a known high-z QSO. The 
quasar J1335-I-3533 can be an optimal candidate for that, 
as it is suitably located close to the SDF. We will consider 
these follow-up ideas in future work. 



APPENDIX A: FULL RADIATIVE TRANSFER 
VS. EXPONENTIAL MODELS 

IZheng et al.] l|2010l ) have studied LAEs using full Lya ra- 
diative transfer in cosmological reionization simulations. 
They compare their results w i th p r evious works (e.g . 
iMcQuinn et al.1 l2007l : llhev et al.l l2008l : iDaval et al] 1201 ll ). 
where the transfer of Lya is treated in an approximate man- 
ner: in such works, in fact, the intrinsic line profile is damped 
by e"'^" , where Ti, i s the scattering optic al depth at frequency 
u along the LOS. IZheng et al.l (|2010l ) emphasize that the 
e"'^" model does not provide an accurate description of the 
observed Lya spectrum, because it does not account for fre- 
quency (no frequency change occurs for any Lya photon) 
and space (no surface brightness information) diffusion. To 
investigate further this point, in addition to the full radia- 
tive transfer model discussed so far, we have reconstructed 

http://www.sdss.org/ 
^ http://sundog.stsci.edu 
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Figure Al. Comparison of the emerging Lya spectra computed 
with the full RT (black) vs. exponential (red) models from the 
central pixel of a simulated prototypical LAE at z = 5.7. The 
dotted line shows the intrinsic line profile. 

the LAE spectra also with the exponential model. We can 
then evaluate the differences between the two. 

The black curve in Fig lAll shows the spectrum from the 
central pixel of the prototypical LAE at z = 5.7 presented 
in Sec. 13.31 The apparent shift of the peak with respect to 
the intrinsic profile is about 5 A. This is consistent with the 
result of Zh eng et al.l (|201Cll ) (see their Fig. 6) for the same 
halo mass, log{Mh) > 10.7 h~^MQ. The red curve in the 
same Figure has been obtained from the e"^" model. 

We find that at all frequencie s such model underesti- 
mates the fiux, but different from I Zheng et al.l (|201(]| ). we 
find the difference between the two methods on 
be < 1%. Most of the scatterings occurring in the cen- 
tral pixel shift the photon frequency in the range within 
8148 A < A < 8153 A, corresponding to the line core and 
thus they are efficiently damped. Average over all the LAEs 
identified in the simulation, the transmissivity computed 
with full radiative transfer is larger by about 5%; some halos 
(0.5% of total) are embedded in a dense HI region, where 
nm ^ 10~'*cm~^. Full radiative transfer of these halos gives 
J.LAE ^ Q g _ while the e~'^" model suppresses the line 
almost completely. 

We have analyzed the possi ble reasons for the discrep- 
ancy between lZheng et al.l (|201G| ) and our findings. We came 
to the conclusion that it originates from the combination of 
two factors: the d ifferent line width and IGM properties. 
IZheng et al.] (|201G| ) used line widths 

^ = 10-*[Mn/{w'°h-'Me)]'^'. (Al) 

This is about 3 times smaller than what adopted here. With 
this line width, almost all the photons are initially within 
the core of the line profile and the fiux of the e"^" model 
is greatly suppressed; as a result only a very small amount 
of photons can escape without any scattering and the main 
contribution to the observed fiux is coming from scattered 
photons. In our case, on the other hand, most photons with 
A > 8155 A are transmitted without any scattering and 
only a very small contribution is coming from scattered pho- 
tons. Even when we enlarged the integrating flux area to 
r ^ 50 kpc the difference with respect to the e"^" model 



remains < 5%. The very small contribution from scattered 
photons is related to a different HI density field. As dis- 
cussed above we have carefully calibrated the IGM trans- 
missivity through an iterative optimization with QSO ab- 
sorption spectra. This proce dure gives on average a 4 times 
lower (riHi) than adopted by IZheng et al.l (|2010l ). which re- 
sults in a reduced spatial/frequency diffusion. 

Finally, the neutral density at the very center of halo 
can be underestimated due to the lack of mass resolution. 
The cell size of L y a rad iative transfer is 1.6 times smaller 
than IZheng et akl (|2010l ) but the mass resolution is 1700 
times larger than theirs. Hence, we have less particles per 
halo, potentially leading to underestimating the recombi- 
nation rate and nm. Tq^^Is determined by the amount of 
transmitted photons rather than scattered photons. How- 
ever, we verify that even with 10 times higher H i density 
in the central cell, the photons with A > 8155 A are trans- 
mitted without scattering since they are far from the core of 
the line profile. The fraction of photons transmitted without 
scattering events is rather sensitive to the peculiar velocity 
of the infalling gas. As we shown in Fig[Sl however, the in- 
fall profile is well captured by remapping SPH particles on 
a finer grid. 
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